Piecing the barcoding puzzle of Palearctic water frogs (Pelophylax) sheds light on amphibian biogeography and global invasions.

Palearctic water frogs (genus Pelophylax) are an outstanding model in ecology and evolution, being widespread, speciose, either threatened or threatening to other species through biological invasions, and capable of siring hybrid offspring that escape the rules of sexual reproduction. Despite half a century of genetic research and hundreds of publications, the diversity, systematics and biogeography of Pelophylax still remain highly confusing, in no small part due to a lack of correspondence between studies. To provide a comprehensive overview, we gathered >13,000 sequences of barcoding genes from >1700 native and introduced localities and built multigene mitochondrial (~17 kb) and nuclear (~10 kb) phylogenies. We mapped all currently recognized taxa and their phylogeographic lineages (>40) to get a grasp on taxonomic issues, cyto-nuclear discordances, the genetic makeup of hybridogenetic hybrids, and the origins of introduced populations. Competing hypotheses for the molecular calibration were evaluated through plausibility tests, implementing a new approach relying on predictions from the anuran speciation continuum. Based on our timetree, we propose a new biogeographic paradigm for the Palearctic since the Paleogene, notably by attributing a prominent role to the dynamics of the Paratethys, a vast paleo-sea that extended over most of Europe. Furthermore, our results show that distinct marsh frog lineages from Eastern Europe, the Balkans, the Near East, and Central Asia (P. ridibundus ssp.) are naturally capable of inducing hybridogenesis with pool frogs (P. lessonae). We identified 14 alien lineages (mostly of P. ridibundus) over ~20 areas of invasions, especially in Western Europe, with genetic signatures disproportionally pointing to the Balkans and Anatolia as the regions of origins, in line with exporting records of the frog leg industry and the stocks of pet sellers. Pelophylax thus emerges as one of the most invasive amphibians worldwide, and deserves much higher conservation concern than currently given by the authorities fighting biological invasions.

Testing Fraternal Birth Order Effects and Antagonistic Effects for Homosexual Men: Power Comparison of Various Methods.

Research on the biological determinants of male homosexual preference has long realized that the older brother effect (FBOE, i.e., a higher fraternal birth rank of homosexuals) and the antagonist effect (AE, i.e., more fertile women have a higher chance of having a homosexual son) can both generate family data where homosexual men have more siblings and more older siblings than heterosexual men. Various statistical approaches were proposed in the recent literature to evaluate whether the action of FBOE or AE could be discriminated from empirical data, by controlling for the other effect. Here, we used simulated data to formally compare all the approaches that we could find in the relevant literature for their ability to reject the null hypothesis in the presence of a specified alternative hypothesis (tests based on regression, Bayesian modeling, or contingency tables). When testing for the FBOE, the relative performance of the different tests was different depending on the specific function generating the older brother effect. Even if no tests were found to always perform better than the others, some tests performed systematically poorly, and some tests displayed a systematic high rate of type-I error. For testing the AE, the relative performance of the tests was generally not changed across all parameter values assayed, providing a clear ranking of the various proposed approaches. Pros and cons for each candidate test are discussed, taking into consideration power and the rate of type-I error but also practicability, the possibility to control for confounding variables, and to consider alternative hypotheses.

Hidden in the sand: Phylogenomics unravel an unexpected evolutionary history for the desert-adapted vipers of the genus Cerastes.

The desert vipers of the genus Cerastes are a small clade of medically important venomous snakes within the family Viperidae. According to published morphological and molecular studies, the group is comprised by four species: two morphologically similar and phylogenetically sister taxa, the African horned viper (Cerastes cerastes) and the Arabian horned viper (Cerastes gasperettii); a more distantly related species, the Saharan sand viper (Cerastes vipera), and the enigmatic Böhme's sand viper (Cerastes boehmei), only known from a single specimen in captivity allegedly captured in Central Tunisia. In this study, we sequenced one mitochondrial marker (COI) as well as genome-wide data (ddRAD sequencing) from 28 and 41 samples, respectively, covering the entire distribution range of the genus to explore the population genomics, phylogenomic relationships and introgression patterns within the genus Cerastes. Additionally, and to provide insights into the mode of diversification of the group, we carried out niche overlap analyses considering climatic and habitat variables. Both nuclear phylogenomic reconstructions and population structure analyses have unveiled an unexpected evolutionary history for the genus Cerastes, which sharply contradicts the morphological similarities and previously published mitochondrial approaches. Cerastes cerastes and C. vipera are recovered as sister taxa whilst C. gasperettii is a sister taxon to the clade formed by these two species. We found a relatively high niche overlap (OI > 0.7) in both climatic and habitat variables between C. cerastes and C. vipera, contradicting a potential scenario of sympatric speciation. These results are in line with the introgression found between the northwestern African populations of C. cerastes and C. vipera. Finally, our genomic data confirms the existence of a lineage of C. cerastes in Arabia. All these results highlight the importance of genome-wide data over few genetic markers to study the evolutionary history of species.

The overlooked importance of vagrancy in ecology and evolution.

Vagrancy is the occurrence of individuals outside the normal geographic range of their species. These rare and unpredictable events have long been neglected by the scientific community, belying a growing body of evidence that vagrancy can have an important role in eco-evolutionary processes at both population and community scales.

Cytogenetics of the Hybridogenetic Frog Pelophylax grafi and Its Parental Species Pelophylax perezi.

Hybrid taxa from the genus Pelophylax can propagate themselves in a modified way of sexual reproduction called hybridogenesis ensuring the formation of clonal gametes containing the genome of only one parental (host) species. Pelophylax grafi from South-Western Europe is a hybrid composed of P. ridibundus and P. perezi genomes and it lives with a host species P. perezi (P-G system). Yet it is unknown, whether non-Mendelian inheritance is fully maintained in such populations. In this study, we characterize P. perezi and P. grafi somatic karyotypes by using comparative genomic hybridization, genomic in situ hybridization, fluorescent in situ hybridization, and actinomycin D-DAPI. Here, we show the homeology of P. perezi and P. grafi somatic karyotypes to other Pelophylax taxa with 2n = 26 and equal contribution of ridibundus and perezi chromosomes in P. grafi which supports F1 hybrid genome constitution as well as a hemiclonal genome inheritance. We show that ridibundus chromosomes have larger regions of interstitial (TTAGGG)n repeats flanking the nucleolus organizing region on chromosome no. 10 and a high quantity of AT pairs in the centromeric regions. In P. perezi, we found species-specific sequences in metaphase chromosomes and marker structures in lampbrush chromosomes. Pericentromeric RrS1 repeat sequence was present in perezi and ridibundus chromosomes, but the blocks were stronger in ridibundus. Various cytogenetic techniques applied to the P-G system provide genome discrimination between ridibundus and perezi chromosomal sets. They could be used in studies of germ-line cells to explain patterns of clonal gametogenesis in P. grafi and broaden the knowledge about reproductive strategies in hybrid animals.

Hybridization in late stages of speciation: Strong but incomplete genome-wide reproductive isolation and 'large Z-effect' in a moving hybrid zone.

In organisms reproducing sexually, speciation occurs when increasing divergence results in pre- or post-zygotic reproductive isolation between lineages. Studies focusing on reproductive isolation origin in early stages of speciation are common and many rely on genomic scans to infer introgression providing limited information on the genomic architecture of reproductive isolation long-term maintenance. This study analyses a natural hybrid zone between two species in a late stage of speciation. We used ddRADseq genotyping in the contact between Podarcis bocagei and P. carbonelli to examine admixture extent, analyse hybrid zone stability and assess genome-wide variation in selection against introgression. We confirmed strong but incomplete reproductive isolation in a bimodal hybrid zone. New findings revealed population genetic structure within P. carbonelli in the contact zone; geographical and genomic clines analysis suggested strong selection against gene flow, but a relatively small proportion of the loci can introgress, mostly within the narrow contact zone. However, geographical clines revealed that a few introgressed loci show signs of potential positive selection, particularly into P. bocagei. Geographical clines also detected a signal of hybrid zone movement towards P. bocagei distribution. Genomic cline analysis revealed heterogeneous patterns of introgression among loci within the syntopy zone, but the majority maintain a strong association with the genomic background of origin. However, incongruences between both cline approaches were found, potentially driven by confounding effects on genomic clines. Last, an important role of the Z chromosome in reproductive isolation is suggested. Importantly, overall patterns of restricted introgression seem to result from numerous strong intrinsic barriers across the genome.

Widespread introgression of MHC genes in Iberian Podarcis lizards.

Major histocompatibility complex (MHC) genes are crucial for the adaptive immune response of jawed vertebrates. Their variation, reaching extreme levels, is driven mainly by an arms race between hosts and pathogens. One hypothesised mechanism contributing to MHC polymorphism is adaptive introgression, the exchange of genetic variants between hybridising species favoured by selection, yet its effect on MHC variation is poorly understood. Detection of adaptive MHC introgression, though challenging, may be facilitated by the analysis of species complexes forming multiple hybrid zones. Here, we investigated MHC introgression in six hybrid zones formed by seven species of Podarcis lizards inhabiting the Iberian Peninsula. To differentiate adaptive introgression from neutral introgression, we compared the patterns of gene exchange in MHC and genome-wide markers. We found elevated sharing of MHC alleles in the proximity of contact beyond the areas of detectable genome-wide admixture in most hybrid zones and, in half of them, asymmetric MHC exchange. In general, the elevated MHC allele sharing between species pairs with abutting ranges compared to geographically isolated species pairs also supports the prevalence of introgression. Collectively, our results demonstrate widespread MHC introgression in the Iberian Podarcis complex and suggest its adaptiveness. Contrary to previous results from Triturus newts, we did not observe differences in the rate of introgression between MHC classes. Our work adds support to the emerging view of adaptive introgression as a key mechanism shaping MHC diversity. It also raises questions about the effect of elevated MHC variation and factors leading to the asymmetry of adaptive introgression.

Genomic phylogeography illuminates deep cyto-nuclear discordances in midwife toads (Alytes).

The advent of genomic methods allows us to revisit the evolutionary history of organismal groups for which robust phylogenies are still lacking, particularly in species complexes that frequently hybridize. In this study, we conduct RAD-sequencing (RAD-seq) analyses of midwife toads (genus Alytes), an iconic group of western Mediterranean amphibians famous for their parental care behavior, but equally infamous for the difficulties to reconstruct their evolutionary history. Through admixture and phylogenetic analyses of thousands of loci, we provide a comprehensive phylogeographic framework for the A. obstetricans complex, as well as a fully resolved phylogeny for the entire genus. As part of this effort, we carefully explore the influence of different sampling schemes and data filtering thresholds on tree reconstruction, showing that several, slightly different, yet robust topologies may be retrieved with small datasets obtained by stringent SNP calling parameters, especially when admixed individuals are included. In contrast, analyses of incomplete but larger datasets converged on the same phylogeny, irrespective of the reconstruction method used or the proportion of missing data. The Alytes tree features three Miocene-diverged clades corresponding to the proposed subgenera Ammoryctis (A. cisternasii), Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), and Alytes (A. obstetricans complex). The latter consists of six evolutionary lineages, grouped into three clades of Pliocene origin, and currently delimited as two species: (1) A. almogavarii almogavarii and A. a. inigoi; (2) A. obstetricans obstetricans and A. o. pertinax; (3) A. o. boscai and an undescribed taxon (A. o. cf. boscai). These results contradict the mitochondrial tree, due to past mitochondrial captures in A. a. almogavarii (central Pyrenees) and A. o. boscai (central Iberia) by A. obstetricans ancestors during the Pleistocene. Patterns of admixture between subspecies appear far more extensive than previously assumed from microsatellites, causing nomenclatural uncertainties, and even underlying the reticulate evolution of one taxon (A. o. pertinax). All Ammoryctis and Baleaphryne species form shallow clades, so their taxonomy should remain stable. Amid the prevalence of cyto-nuclear discordance among terrestrial vertebrates and the usual lack of resolution of conventional nuclear markers, our study advocates for phylogeography based on next-generation sequencing, but also encourages properly exploring parameter space and sampling schemes when building and analyzing genomic datasets.

The Yellow-browed Warbler (Phylloscopus inornatus) as a model to understand vagrancy and its potential for the evolution of new migration routes.

Why and how new migration routes emerge remain fundamental questions in ecology, particularly in the context of current global changes. In its early stages, when few individuals are involved, the evolution of new migration routes can be easily confused with vagrancy, i.e. the occurrence of individuals outside their regular breeding, non-breeding or migratory distribution ranges. Yet, vagrancy can in theory generate new migration routes if vagrants survive, return to their breeding grounds and transfer their new migration route to their offspring, thus increasing a new migratory phenotype in the population. Here, we review the conceptual framework and empirical challenges of distinguishing regular migration from vagrancy in small obligate migratory passerines and explain how this can inform our understanding of migration evolution. For this purpose, we use the Yellow-browed Warbler (Phylloscopus inornatus) as a case study. This Siberian species normally winters in southern Asia and its recent increase in occurrence in Western Europe has become a prominent evolutionary puzzle. We first review and discuss available evidence suggesting that the species is still mostly a vagrant in Western Europe but might be establishing a new migration route initiated by vagrants. We then list possible empirical approaches to check if some individuals really undertake regular migratory movements between Western Europe and Siberia, which would make this species an ideal model for studying the links between vagrancy and the emergence of new migratory routes.

DNA barcode reference library for the West Sahara-Sahel reptiles.

DNA barcode reference libraries are now continuously produced for the tree of life, which are essential pillars for the study of biological diversity. Yet, our knowledge about global diversity is largely limited in undersampled regions such as the largest warm desert, the Sahara-Sahel. This dataset provides a DNA barcode reference library for the reptiles of the Western Sahara-Sahel (WSS) and neighbouring countries across this region. It includes 760 barcodes from 133 reptile taxa, distributed in 23 families, and covering the intraspecific diversity of some species. A total of 84 species were collected in the WSS (83% of the total reptile species richness) over 18 overland field expeditions conducted since 2003. DNA barcodes resulted in a high success rate (95%) of species identification and barcoding gap analysis highlighted the effectiveness of the COI fragment as a barcode marker for the WSS reptiles. This dataset represents a comprehensive and reliable DNA reference library for the WSS, filling an important biodiversity gap across a remote and hard-to-sample region.

Sex chromosomes as supergenes of speciation: why amphibians defy the rules?

As reflected by the two rules of speciation (Haldane's rule and the large X-/Z-effect), sex chromosomes are expected to behave like supergenes of speciation: they recombine only in one sex (XX females or ZZ males), supposedly recruit sexually antagonistic genes and evolve faster than autosomes, which can all contribute to pre-zygotic and post-zygotic isolation. While this has been mainly studied in organisms with conserved sex-determining systems and highly differentiated (heteromorphic) sex chromosomes like mammals, birds and some insects, these expectations are less clear in organismal groups where sex chromosomes repeatedly change and remain mostly homomorphic, like amphibians. In this article, we review the proposed roles of sex-linked genes in isolating nascent lineages throughout the speciation continuum and discuss their support in amphibians given current knowledge of sex chromosome evolution and speciation modes. Given their frequent recombination and lack of differentiation, we argue that amphibian sex chromosomes are not expected to become supergenes of speciation, which is reflected by the rarity of empirical studies consistent with a 'large sex chromosome effect' in frogs and toads. The diversity of sex chromosome systems in amphibians has a high potential to disentangle the evolutionary mechanisms responsible for the emergence of sex-linked speciation genes in other organisms. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

A lot of convergence, a bit of divergence: Environment and interspecific interactions shape body colour patterns in Lissotriton newts.

Coexistence with related species poses evolutionary challenges to which populations may react in diverse ways. When exposed to similar environments, sympatric populations of two species may adopt similar phenotypic trait values. However, selection may also favour trait divergence as a way to reduce competition for resources or mates. The characteristics of external body parts, such as coloration and external morphology, are involved to varying degrees in intraspecific signalling as well as in the adaptation to the environment and consequently may be diversely affected by interspecific interactions in sympatry. Here, we studied the effect of sympatry on various colour and morphological traits in males and females of two related newt species Lissotriton helveticus and L. vulgaris. Importantly, we did not only estimate how raw trait differences between species respond to sympatry, but also the marginal responses after controlling for environmental variation. We found that dorsal and caudal coloration converged in sympatry, likely reflecting their role in adaptation to local environments, especially concealment from predators. In contrast, aspects of male and female ventral coloration, which harbours sexual signals in both species, diverged in sympatry. This divergence may reduce opportunities for interspecific sexual interactions and the associated loss of energy, suggesting reproductive character displacement (RCD). Our study emphasizes the contrasting patterns of traits involved in different functions and calls for the need to consider this diversity in evolutionary studies.

A new westward migration route in an Asian passerine bird.

The evolution of migration routes in birds remains poorly understood as changes in migration strategies are rarely observed on contemporary timescales.1-3 The Richard's Pipit Anthus richardi, a migratory songbird breeding in Siberian grasslands and wintering in Southeast Asia, has only recently become a regular autumn and winter visitor to western Europe. Here, we examine whether this change in occurrence merely reflects an increase in the number of vagrants, that is, "lost" individuals that likely do not manage to return to their breeding grounds, or represents a new migratory strategy.4-6 We show that Richard's Pipits in southwestern Europe are true migrants: the same marked individuals return to southern France in subsequent winters and geo-localization tracking revealed that they originate from the western edge of the known breeding range. They make an astonishing 6,000 km journey from Central Asia across Eurasia, a very unusual longitudinal westward route among Siberian migratory birds.7,8 Climatic niche modeling using citizen-science bird data suggests that the winter niche suitability has increased in southwestern Europe, which may have led to increased winter survival and eventual successful return journey and reproduction of individuals that initially reached Europe as autumn vagrants. This illustrates that vagrancy may have an underestimated role in the emergence of new migratory routes and adaptation to global change in migratory birds.9,10 Whatever the underlying drivers and mechanisms, it constitutes one of the few documented contemporary changes in migration route, and the first longitudinal shift, in a long-distance migratory bird.

Mass of genes rather than master genes underlie the genomic architecture of amphibian speciation.

The genetic architecture of speciation, i.e., how intrinsic genomic incompatibilities promote reproductive isolation (RI) between diverging lineages, is one of the best-kept secrets of evolution. To directly assess whether incompatibilities arise in a limited set of large-effect speciation genes, or in a multitude of loci, we examined the geographic and genomic landscapes of introgression across the hybrid zones of 41 pairs of frog and toad lineages in the Western Palearctic region. As the divergence between lineages increases, phylogeographic transitions progressively become narrower, and larger parts of the genome resist introgression. This suggests that anuran speciation proceeds through a gradual accumulation of multiple barrier loci scattered across the genome, which ultimately deplete hybrid fitness by intrinsic postzygotic isolation, with behavioral isolation being achieved only at later stages. Moreover, these loci were disproportionately sex linked in one group (Hyla) but not in others (Rana and Bufotes), implying that large X-effects are not necessarily a rule of speciation with undifferentiated sex chromosomes. The highly polygenic nature of RI and the lack of hemizygous X/Z chromosomes could explain why the speciation clock ticks slower in amphibians compared to other vertebrates. The clock-like dynamics of speciation combined with the analytical focus on hybrid zones offer perspectives for more standardized practices of species delimitation.

Nuclear phylogenies and genomics of a contact zone establish the species rank of Podarcis lusitanicus (Squamata, Lacertidae).

Unravelling when divergent lineages constitute distinct species can be challenging, particularly in complex scenarios combining cryptic diversity and phylogenetic discordances between different types of molecular markers. Combining a phylogenetic approach with the study of contact zones can help to overcome such difficulties. The Podarcis hispanicus species complex has proven to be prosperous in independent evolutionary units, sometimes associated with cryptic diversity. Previous studies have revealed that one of the species of this complex, P. guadarramae, comprises two deeply divergent yet morphologically indistinguishable evolutionary units, currently regarded as subspecies (P. g. guadarramae and P. g. lusitanicus). In this study we used molecular data to address the systematics of the two lineages of Podarcis guadarramae and the closely related P. bocagei. Firstly, we reconstructed the species tree of these three and two additional taxa based on 30 nuclear loci using the multispecies coalescent with and without gene flow. Secondly, we used SNPs obtained from RADseq data to analyze the population structure across the distribution limits P. g. lusitanicus and P. g. guadarramae, and for comparison, a contact zone between P. bocagei and P. g. lusitanicus. Nuclear phylogenetic relationships between these three taxa are clearly difficult to determine due to the influence of gene flow, but our results give little support to the monophyly of P. guadarramae, potentially due to a nearly simultaneous divergence between them. Genetic structure and geographic cline analysis revealed that the two lineages of P. guadarramae replace each other abruptly across the sampled region and that gene flow is geographically restricted, implying the existence of strong reproductive isolation. Podarcis bocagei and P. g. lusitanicus show a similar degree of genetic differentiation and reproductive isolation, with very low levels of admixture in syntopy. These results support that all three forms are equally differentiated and reproductively isolated. In consequence, we conclude that the two former subspecies of Podarcis guadarramae constitute valid, yet cryptic species, that should be referred to as P. lusitanicus and P. guadarramae.

Systematics, biogeography and evolution of the Saharo-Arabian naked-toed geckos genus Tropiocolotes.

Plate tectonics constitute one of the main mechanisms of biological diversification on Earth, often being associated with cladogenetic events at different phylogenetic levels, as well as with exchange of faunas and floras across previously isolated biogeographic regions. North Africa and Arabia share a complex geological history that dates back to the break-up of the Arabian plate from the African plate ~30-25 Mya, followed by various geological events, such as the formation of the Red Sea or the connection between the African, Arabian and Eurasian plates. Species with Saharo-Arabian distributions have shown a close association between their evolutionary history and these geological events. In this study, we investigate the systematics, biogeography and evolution of the genus Tropiocolotes, a group of small ground-dwelling geckos, comprised by 12 species distributed from the Atlantic coast of North Africa to southwestern Iran. Species delimitation analyses uncovered the existence of high levels of undescribed diversity, with forms here considered at the species level including Tropiocolotes tripolitanus (Mauritania and southern Morocco), T. nattereri (southern Israel) and T. scorteccii (Yemen and Oman). Phylogenetic and biogeographic analyses recovered two main clades, an exclusively African clade and a Saharo-Arabian clade, that split ~25 Mya following the vicariant event mediated by the separation of the Arabian and African plates. The complex geological activity around the Red Sea is associated with the diversification within the Saharo-Arabian clade, including the colonization of North Africa from a second Tropiocolotes group. Results also provide new insights into the geographic distribution of Tropiocolotes nubicus, previously considered as exclusively associated to the Nile River valley, extending its known distribution further west, up to the Central Mountains of the Sahara. Accordingly, the Nile River seems to act as a major biogeographic barrier, separating Tropiocolotes nubicus and T. steudneri in their western and eastern margins, respectively.

Variable levels of introgression between the endangered Podarcis carbonelli and highly divergent congeneric species.

Recent empirical studies have demonstrated that speciation with gene flow is more common than previously thought. From a conservation perspective, the potential negative effects of hybridization raise concerns on the genetic integrity of endangered species. However, introgressive hybridization has also been growingly recognized as a source of diversity and new advantageous alleles. Carbonell's wall lizard (Podarcis carbonelli) is an endangered species whose distribution overlaps with four other congeneric species. Our goal here was to determine whether P. carbonelli is completely reproductively isolated from its congeners and to evaluate the relevance of hybridization and interspecific gene flow for developing a conservation plan. We used restriction site associated DNA (RAD) sequencing to discover SNPs in samples from four contact zones between P. carbonelli and four other species. Principal component analysis, multilocus genotype assignment and interspecific heterozygosity suggest incomplete reproductive isolation and ongoing gene flow between species. However, hybridization dynamics vary across all pairs, suggesting complex interactions between multiple intrinsic and extrinsic barriers. Despite seemingly ubiquitous interspecific gene flow, we found evidence of strong reproductive isolation across most contact zones. Instead, indirect effects of hybridization like waste of reproductive effort in small isolated populations may be more problematic. Our results highlight the need to further evaluate the consequences of introgression for P. carbonelli, both on a geographic and genomic level and included in a comprehensive and urgently needed conservation plan. Besides, those findings will add important insights on the potential effects of hybridization and introgression for endangered species.